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Ansj'KAn; 

Methods  of  modelling  probabi 1 is t ie  components 
which  are  not  adequately  represented  hv  the  .stan- 
dard continuous  distributions  (such  as  norn.il,  gam- 
ma, and  Wei  bull)  art’  surveyed.  The  methods  are 
categorized  as  systems  of  distributions,  approxi— 
mations  to  the  cumulative  distribution  (unction, 
and  four-parameter  distributions.,.  . basis  is  on 
generality,  determin.it  ion  •:  appropriate  par.itr.eter 
values,  and  random  variate  generation. 


Johnson 

Approximations  of  the  inverse  cdf 
Kxp.ins ion  techniques 
Polynomial  regress  ion 
Kectangu 1 ar  approximati on 
Four-parameter  d i st  ri but  ions 
Beta 

Fo  u r - p a rank*  ter  gamma 
Burr 

Generalized  lambda 
Absolute  lambda 


±. INTUODICTION 

A common  problem  in  computer  simulation  of  stochas- 
tic systems  is  t..-  selection  of  appropriate  statis- 
tical distributions  tor  ,*:n  rating  random  observn- 
( • fh<  elect i • . • in  over-si  H- 

f ied  process  of  s elect  in;:  from  a few  well-known 
standard  dJ  tri  it  ion  : i ri  il , to  • n il , - nuna, 

and  Wei bull  (including  the  exponential)  in  the  con- 
tinuous s . i ; Bernoulli,  binomial,  and  Poisson  in 
the  discrete  case.  These  distributions  are  appro- 
priately chosen  in  many  cases.  Other  times  they 
arc  chosen  because  of  the  ease  of  their  pseudo- 
random process  generation. 

Ideally,  an  appropriate  statistical  distribution  is 
chosen  using  a three  step  iterative  process  consis- 
ting of 

1.  identification  of  a family  of  distribu- 
tions which  appear  appropriate 

2.  determination  of  the  parameter  values 
which  allow  the  family  to  best  fit  the 
a itu.it  ion,  and 

3.  diagnostic  checking  to  determini  whether 
or  not  an  adequate  model  has  been  found. 

These  three  steps  are  repeated  until  the  diagnostic 
check  is  satisfactory.  (McGrath  et  al  [28 j give 
a detailed  survey  of  techniques  for  selecting  pro- 
bab i 1 i t y dic.tr  i but  ions. ) 

In  those  cases  where  the  diagnostic  cheek  is  not 
satisfactory  for  the  standard  distribution,  more 
elaborate  means  must  be  u ied.  This  paper  surveys 
methods  useful  r»r  continuous  random  variables. 

In  particular  the  following  are  discussed: 

Systems  of  Distributions 

Pearson 


The  common  property  of  these  methods  is  their  abil- 
ity to  satisfy  four  or  more  desired  proport ie*  . 

The  standard  distributions  usually  satisfy  one 
(exponential),  two  (normal,  gamma,  Weibull),  «*r 
three  (translated  gamma  or  Weibull)  properties. 

These  proper  t ie.*;  may  be  moments,  quantiles,  loca- 
tion of  the  mode  ‘or  other  desired  characteristics. 
The  generality  of  n method  depends  on  the  number  of 
properties  satisfied  and  the  combinations  of  pro- 
pert  i es  pos s ib 1 e . 

The  generality  of  an  approach  nay  be  measured  to 
some  extent  by  considering  the  third  and  fourth 
standard! zed  moments 

aj  = K{  [ (X-O/o  ] 1 ) i = 3 . 4 . 

where  u is^lhe  mean  and  ' the  standard  deviation  of 
X.  3.  ~ Ot  is  a measure  of  skewness  and  the  kui  t - 
sis  (tail-weight  and/or  peakodness)  may  be  measured 
by  8^  = ot^ . Since  (p  , p0)  is  a measure  of  the 
shape  of  a distribution,  the  more  values  of  . ,) 

attainable  by  a method,  the  more  general  it  is. 

While  other  measures  are  also  reasonable,  standar- 
dized moments  will  be  used  here  t e motivate  these 
methods  and  to  compare  their  generality. 

Figure  1 shows  the  (•  j»  8,)  plane.  The  locutions 
of  Bernoulli  trial-.,  uniform,  normal,  Weibull,  gam- 
ma, exponential  and  St udent-t  are  shown.  Note  that 
symmetric  d i tribulions  lie  on  the  8.  0 axis. 

Bernoulli  trials,  the  shortest  tailed  distribution 
possible,  lie  oti  the  ~ B j 1 line.  The  ?.\mr*e- 
tric  Bernoulli  trial  (p  = . '>)  lies  at  the  point 
( 8 j . 8 , ) (0* 1) . i he  norn jl  distribution  lie  at 

(0,1)  and  Student's  t distribution  lies  below  the 
normal.  Asymmetric  distributions  are  scat tt ied 
around  the  square.  The  exponential  is  at  (:  , ; ,) 

K (4,9)  which  is  a point  on  both  the  Weihull  and* 
gamma  curves.  Hahn  and  Shapiro  | 1 b | discus*,  tin 
positions  of  the  standard  distributions  in  r.ote 
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A characteristic  of  all  of  those  distributions  is 
that  none  covers  a revicn  of  the  (B  , B.,)  plane. 

The  dist  r ibut  ions  usually  considered  general  (Wei- 
bull,  lognormal,  gamma  and  Student's  t)  are  asso- 
ciated with  lines  while  more  limited  models  (nor- 
mal, exponential,  unil\> rm)  are  represented  by 
points.  It  is  clear  that  the  standard  distribu- 
tions correspond  to  only  a small  part  of  the  .avail- 
able models  of  randomness.  While  the  standard 
techniques  are  often  appropriate  din  to  their  theo- 
retical foundation,  when  theoretical  reasons  do  not 
justify  a standard  distribution,  selection  f rom 
any  point  on  the  plane  is  reasonable  and  often  de- 
sirable. 

The  property  common  to  all  the  versatile  methods 
mentioned  at  the  beginning  of  this  see t ion  is  that 
each  covers  a region  of  the  (B. , B0)  plane.  To 
have  thi-.  versatility,  at  I east  four  parameters  are 
necessary:  two  for  shape  (third  and  fourth  moments), 

one  for  scaling  (second  moment  --  variance)  and  one 
for  location  (first  moment  — mean).  Some  models 
have  more  parameters. 

The  three  groups  of  methods  are  examined  in  the 
next  three  sections.  Generality,  properties,  tech- 
nique «»t  parameter  determination  and  techniques  of 
random  variate  generation  are  discussed.  The  au- 
thor apologizes  for  any  omissions  and  would  appre- 
ciate having  them  brought  to  hi.,  attention. 

i i . S-,  sn.MS  or  ntsi'Kip.i  iions 

The  Pearson  and  .Joitn  -on  systems  each  cover  the  en- 
tile (*  j»  •' ,)  plane  l>v  using  multiple  functional 
forms,  bach  functional  form  covers  a itfiim  * * f the 


(B. , B ,)  plane  or  is  a transition  between  regions. 
Both  sy-ter.s  were  developed  as  statistical  models 
before  computer  simulation  was  a factor. 

PEARSON  SYSTEM 

The  Pearson  system  (321  is  composed  of  seven  types 
(functional  forms).  Some  correspond  to  standard 
distributions  (beta,  gamma.  Student  t)  and  the 
others  provide  a continuum  over  the  (f*  , B^)  plane. 
There  is  exactly  one  Pear-. on  distribution  Tor  each 
point  of  the  (i-  j , B?)  plane.  The  normal  distribu- 
tion is  a limiting  <!  i s t r i but  i on  of  every  type  of 
Pearson  distribution.  Elderton  and  Johnson  |13], 
Johnson  and  Kotz  [24]  and  Ord  [30]  discus-,  the 
Pearson  system  in  detail.  Tadikanalla  | 36 ] dis- 
cusses type  selection,  parameter  estimation,  and 
^iriate  generation. 

Due  to  the  multiple  forms,  multiple  parameter  esti- 
mation and  variate  generation  techniques  are  neces- 
sary to  cover  the  whole  ( .• , , r\,)  plane.  Johnson, 
et  nl  [23]  give  tables  of  percentage  points  of  the 
Pearson  curves  for  given  ».  . and  r , which  can  aid 
parame ter  determination.  Max in 

mat  ion  is  difficult.  However,  parameters  may  be 
determined  by  matching  the  first  four  moments  rela- 
tively easily.  The  Pearson  system  is  on«-  ef  the 
few  approaches  having  closed  form  solution  tor  mat- 
ching moments.  For  example,  Pearson  tv;vs  I and  II 
(beta  distribution)  yield  parameter  value-. 

P,q  = r{ l+(r+2)  i&  {(r+2)  ? +lf.(r+l))  ) }/2 

where 

r = 6(R2-61-l)/(6+361-2f:2) 

and  p < q according  to  whether  /3  g 0.  The  onlv 
drawback  is  that  the  equations  change  tor  each  type. 
However,  the  type  may  he  found  easily  from  ( b ^ , . .) 
charts  usually  reproduced  whenever  the  Pearson  sy- 
stem is  discussed. 

Variate  generation  is  not  straightforward.  Each 
type  requires  its  own  techniques.  The  inverse  cdf 
is  not  closed  form  for  the  more  important  types  and 
rejection  techniques  have  not  worked  well  until 
recently  due  to  many  types  having  long  tail;..  (cop- 
er, Davis,  and  Done  ( 0 ] do:.,  ribo  theii  1TK(T  pro- 
gram which  calculates  moments  and  generates  vale's 
using  numerical  integration.  Numerical  integral  ion 
was  important  at  that  tire  since  no  exact  techni- 
ques were  known.  Since  then  much  development  has 
occurred  in  generation  techniques  Jor  the  IV  it  . on 
system.  In  particular,  tvp.  I and  11  Peat  - on  (beta 
variate  general  ion  lias  been  the  subject  tor  Johnk 
1 19 1 , Fishman  | l > 1 , Whittaker  (All,  Ahrens  and 
Dieter  (1),  Atkinson  and  P»  arce  ( | , Atkinson  and 
Whittaker  | 1],  and  S ltmeiser  and  Shalahv  | 5f>|. 

These  developments  are  important  to  the  whole  Pear- 
son family  since  beta  vat  i tics  can  he  transformed 
to  obtain  either  types  difficult  to  generate  direct- 
ly. Psing  the  beta  ro lat ionship , V ‘-rath  and  Ir- 
ving, [2‘1]  give  flow  charts  and  FORTRAN  subroutines 
for  each  Pearson  tvpo. 
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John  . on  I/O,  J 1 1 has  dew  1 . ;u*«i  two  sv-.tir  . which 
cow i t In*  • lit  it*'  ( ; | » ; , ) p l .im-  in  tin..  1 mu  t i on- 
al  torus  bused  on  t t an  ■ 1 o t it  lens  ot  tin  noim.il 
and  I ipl.iiu  Ji-.tribut  ion...  l ot  .•  be  a st und.ird 
normal  v.ir lute.  Thin 


x a (,  t exp l 

r k ( : to  exp  | ( 'I  1 1 
X It  exp 1 1 z-  , I 1 


1 ’l  < • ) • l ■ I 


y - Zj  - in  ( x t - r ) 

yield:,  two  di  .ivril  percent  ill 
al  so  d iscus-  es  t lu  -e  t op  its. 


Tad  i karna  1 1 a | 0*  ] 
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x *=  t.  + 1 si nli  l (z.~y)/  ‘ J < x • •• 

correspond  to  tin  three  types.  S(  ( 1 ognorma l ) , S 
(hounded)  and  S^.  (unbounded),  yielding,  straight 
forward,  though  relatively  slow,  variate  general  ion. 
The  second  system  is  obtained  bv  letting  .•  be  a 
standard  Laplace  random  variate.  For  both  sv steins 
tr  is  a location  parani*ter,  ■ is  a scale  parameter, 
and  Y and  S are  shape  parameters. 

Parameter  estimation  can  proceed  in  at  least  two 
ways  : mat  chi  tip,  percent  i K . or  matching  moments. 

Johnson  [//]  gives  tables  for  i and  ‘ in  terms  of 
and  ; for  the  S family.  Johnson  and  Kotz 
1 2 * ) giw'uri  abac  with  the  sank*  inl4>rmat  ion  in 
graphical  form.  (Liven  Y and  , the  desired  vari- 
ance 0‘  is  obtained  from 

X * o/o(y) 

and  the  desired  m-*an  p is  obtained  from 
f.  ~ p - X K(y) 

where 

1 / 2 

K(y)  -u)  ' si  nli  Q 


0(y)  l(w-l)(u)  cosh  2iM-l )/  / 1 

-?  ( 
u)  = exp  (6  ‘ ) 


1 1 1 1 . «•  teel  liqtii  basr  i i it  ing  the  Invei 

eumulatiw  distribution  turn  t ion  (cdM  have  b«  on 

proposed s ex  t<  , polyi  lal  regres 

• : ion,  and  roc  t nngu  1 a i approx  it  it..:,.  Al  thou 
three  only  approximate  tin  desired  cdf,  tin-  appie- 
ximation  can  be  made  as  accurate  a*,  desired  at  the 
expense  of  coding,  storage  and  -a  i xecut ive  t < . 


EXPANSION  TKUlNlhUTS 

Expans  ion  techniques  can  sometime;  provide  a random 
variable  having  tout  or  non  disired  rmont  . Ihe 
general  technique  is  1 1’  expand  t:-  den-  itv  or  cdt 
into  a series  of  inf  in  it.-  terms  based  on  some  ini- 
tial distribution  which,  ii  used  in  its  entiretv, 
would  provide  the  desired  moments  exact lv.  Tin-  ini- 
tial distribution  i ••  usually  the  normal,  a 1 1 In  mcIi 
others  siiili  a.-,  tin  gamma  hav«*  been  used.  Eldert on 
and  Johnson  (1  1),  dohns on  and  Kot  .•  I * 1 and  Ord  | 10) 
discuss  expansion  technique:,  in  more  detail.  Only 
the  special  case  of  the  Corn  i sh-  1*  i slier  expansion 
[101  i s di  sells  ed  here  . 

The  Corni sh-F i slier  expansion,  which  uses  a normal 
initial  distribution,  cun  be  used  to  generate  ran- 
dom variates  having  desired  moments  J (mean),  '..,th 
(variance),  (third  central  moment),  ...,  (k 
central  moment).  We  consider  k A here.  Set 


12  « y/6 

For  tin*  S family,  \ is  the  lower  bound  and  x is 
tin*  range,  re<lm*ing.  the  problem  to  only  two  parame- 
ters y and  '.  Oiven  any  two  desired  percentiles 
(x . , Pj ) and  (x,,  p,)  yields 


= p4  " p;> 

which  are  the  first  four  cumul.tnfs,  upon  whicli  the 
lorn ish-Fi slut  expansion  is  based.  Now  generate  j 
normal  (0,1)  variate  Setting 


1 •< 

V " 

]*”i 

x . 

l n 

" N 

1 III 

1 1 

^ f V - x 

a “ lS  111 

V . 
■>x- 

c 

x'j 

/j  ? !(p.),  the  p.  percentile  of  the  stan- 

dard nor n.i  l u i st  r i but  ion . 

Similarly  f * » r the  S distribution,  is  the  lower 


z.  -1  (/.*'- 1 ) ».  j/h  * (.*  - L ■)c4//4 
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yields  v » I , \ as  the  desired  variate.  Fish 
el  and  lolllilli  [1*1  ive  tel:  thlough  the  el  hth 


m i : • • ’ : * • ; v’-;  v. i : is.;  ay*  « \ v • * : i \. 


(‘..lit  i mu'il 


cumul ant . 

Kxpansion  U i ha  iqui  •.  cotwcrp.  rapidly  when  the  dis- 
tribution desired  is  imilat  to  tin*  initial  distri- 
bution, but  later  t erms  luiorn-  important  otherwise, 
(barton  and  ivutu*.  [S]  discus;  tin*  rep  ion  ot  (,  , 

Iv,)  where  the  i."  nit  iii;.'.  distribution  is  unimodal 
and  ha*;  .ill  nonnep.it  ive  din*  ilv.)  Iheietotc  euro 
must  bt*  taken  win  n usinp,  tilt*,  technique.  A simple 
pro.  i dui  i'  is  to  ioll.it  -ample  moments  on  tin*  » n- 
t ratod  ob-.orvat  iisn.  to  onsuro  tin*  desired  moments 
art*  obtained. 


tin  ntinlur  of  If  neat  •«>:.nfs  prow  , and  flier.  j..|. 
tlu*  » rior  dci  riMsis,  tin  dilti  ult%  in  ohtaininp 
de-.irid  mements  i in  i « • . i *; . • *. . Ni-vrrt  ho  loss  , t b i r « < » i — 
niqne  works  rvrn  wlirn  tin-  d.  si  nd  distribution  i- 
not  unimodal,  which  i not  trim  ot  t lie  other  r .»■  t hud*, 
dis.u-.oil  ho  re  . Hut  lor  [ | extends  tret  angular  ap- 

prox in. it  i on  to  tin  equivalent  ot  quadratic  appr.xf- 
Bu»t  ion  by  usinp.  i.jo.  ti<>n  methods  on  a piccevi-e 
1 i him  t di'iis.  i t v t mn  t i >n  . 
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Tocher  1.181  uses  an  approximat  i>':i  to  t lu*  invorso 
oil 

.-1  2 2 2 

K (u)  a+buhu  ti(l-u)  I op.  I u)  ■* . u*  I op,  (,  1 -u ) 
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KI  ( l ANCn.AK  AITKOXINA  IT  ON’ 

Hoc  t anp.ul  ar  approxiri.it  ion  i s pn-babl\  f ho  most  vide- 
lv  u*. od  tcclmiqiu  tor  p.enerat  inp  random  variates 
which  d jot  t f i ■ th<  standard  distri l-u i ion  • 

( , for  es  ■ . us«  -.:!ii  • • hn  ue . • 1 in  i n 

Vi*rso)  I’d  t i--  appt  o\ i ; iat  o.l  bv  a p i oe vw  i • linear 
filin' t i • » n . wliieli  i • equivalent  to  a mixture  ot  k 
nonover lappi np  adjacent  uni  torn  d i • t i i but  i on*. . To 
> ' S ; It  . : ' • . ' ' ' tie  1 

crated  from  tin  i di  tribntiou  with  probability 

I’.,  i 1 , k.  In  the  -peel  il  ease  ot  p 

J/k  lor  all  i , ■ f o r a p e is  redm  «•<!  in.,  tin-  p.’s  no 
Jonp.or  nood  to  In  tor.  i.  Met'.tath  and  ltvinpl|uM 
and  K.irn.itd  and  i awd.*t  v |.|  point  out  that  tin-  con- 
stant p value-;  o.i*.«*  ..in  also  I*,  made  la  . ter  than 
p.eno/al  r oe ; an  ;u  1 a i approximation  when  k i ••  larpe. 

While  variate  penerat Ion  i.  diic.t,  the  proper 
piecewise  lino  it  appt ox imat ion  is  not  so  obvious, 
brake  and  t.Vrhardt  jlll  ami  Met. i lehr i . f |.'/|  s ttidy 
this  problem  in  a p.netal  nontext.  Ktsro  | .'(*  j if  i s- 
cuksm  thi.  problem  in  t ho  simulat  ion  context  * pre- 
sent inp  a 2 7 step  a l cor  1 1 1:*,  to  tit  .sipirieal  data. 
Even  with  an  a ! cor  i f Irrt  i approueh  if  i . dfltieult 
to  oht  liu  t ho  do*;  j red  r.onu  nt  . l ot  . sample,  in 
(Cisco's  exa  ' I I t I ' ' • -oh,  t lie 

me. in  i*.  .Si.'  uni.  • too  lat  *.  . Vhis  m in  ot  ni.tv  not 
be  important  itt  a p i von  situation.  (M  .out  .e  a 
non  I i ti  n pro  t ■ inp  problem  could  1 *.olvod  alonp, 
the  1 in.--,  ol  : t Hi  . i . i n ave  i a , t lot  sub  je.T  to 
sat  i I - il  .,  the  i!e  itril  ti.  an  and  vatiaaeo.  Hut  a*. 


There  ate  five  distributions  which  covet  wide  re- 
gions ot  tin*  U'j.  I ,)  piano:  beta  (altsadv  di  - 

cussed  under  the  lY.'ir  .on  system),  font'  param*  t<  t 
p.'tmma,  lUtrr,  penera 1 i ccd  lambda  d i s t r i but  i on  , and 
a distribution  referred  to  Jjere  as  the  absolute 
lambda  distribution.  All  have  exactly  tout  parame- 
ters and  only  one  luin-t  ii'ii.d  form,  nsnltinp  in 
parameter  estimation  and  variate  general  ion  methods 
which  arc  applicable  in  many  mode  Mini;  situations 
white  the  standard  distributions  are  not  adequate. 


KOUtt- I’AkAMKTKU  tIAMMA  1*  I STKI  Ht’TlON 

Harter  |17J  discusses  a four-pararu-t  er  distribution 
which  includes  several  standard  d i s t r i lui  t i oils  as 
jopceial  cases.  The  density  function  is 

f ( X ; e , a , b , p ) *p(x--c)^^  * exp{- | t,  x -e  ) /a  1 1 : . a' i (b) 

where  a,b,p  N 0 and  x s e.  Here  c is  tin*  location 
parameter,  a is  the  scale  parameter  and  b ami  p : 
pet  her  determine  shape.  When  e 0,  this  distribu- 
tion redueis  to  the  e.nmma  il  p 1,  Weibull  ill*  l, 

exponential  i l b p i,  mid  half-normal  i f b 
pT'.  In  additi<>n,  the  normal  distribution  is  a li- 
miting ease  as  p 1 and  b p.oes  t>*  intinitv. 

While  beiny*  more  perietal  than  either  the  Weibull  or 
pantnt.i  d i s t r Unit  ion , this  distribution  has  ov.i,.] 
d i s.nlv. ml  apes  l ot  use  in  simulation.  First,  the 
moments  are  not  e.is.v  to  evalnite  except  tor  tin* 
special  cases.  Second,  parameter  estimation  i not 
s t r.t  ip.ht  t orward . Harter  (17j  diseussos  maximum 
likelihood  estimation  ineludinp,  the  est  in.itoi  s' 
asvmptotie  v.iriune.s  ami  e.’V.ir  i.inees . The  third 
di  s.uivant  ape  i *;  in  random  variate  p.enerat  ion.  Since 
the  standard  p.amma  distributions  (Teat  son  type  111) 
are  not  s t rai  p.ht  t orward  to  ponerute  fexeept  for  in- 
teger shape  pa  ratno  t «r ) , eortainly  this  tour  pat.  « 
ter  pi  net  .i  1 i /at  ion  i .*•  dillieult  to  I'.i’tU'  rate.  The 
author  knows  ol  no  results  in  thi?;  area. 


HTKJv  PIS  I KlIUTTi^N 


I'lie  Hurt  il  i s t t i but  ion  has  edt 


!(x)  l - | I + {(x  a: 

with  pat ame lets  h and  k po- 
lo, at  ion  t»ar.tmet  er  , b ( he  : 
an.l  k determine  the  shape, 
by  Hut t l s 1 , t he  paramo t er 
shown  bv  Win  elei  [ a'l  the  i 
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by  this  ili-.tribut  in  i tael  mb ■*.  the  W»  tbull  .uul  >• . un— 
m.i  ourws  anil  • at  eu  eon . ••  1 > a ! i • . . !>'  In  av  i » t 
tail--  (l.irget  ,).  Kurr  | ‘1  r.oti  tl.it  u-ini'  . 0 

results  in  c o v » : I » . ' i a now  t * 1*  ::  ! i : j • • , ) 

closet  to  the  line  . . - . j I in  t1-  beta  region. 

Whoilet  1*0]  and  . . ' 1 k . I 1 .1  .nil  !‘  . : • » • t g ( i / ] * , . 1 \t  • 

URt'd  the  Sun  . • 1 . ( . 1 t imal  t 1 

d i st  r i but  ion . 

The  j’.ii  1 etots  mav  bo  estimated  usim:  tables  >•  i \ . 1 

in  Burt  1 7 1 to  nut  eh  the  tit  t tout  '..wnts.  Max- 
imum likelihood  i t r otors  t in  be  obtained  via 
iterative  tm  et  t . a I !y  t echn  1 .,m-s  . 1..  kow  1 • « ttving 

t o tit  the  pa  1 a:\eter  ; bv  spre  i ! s i :.g  quant  i !os  n 
quires  mine!  ical  olut  ion.  Random  vat  sate  genera- 
t ion  is  direel  I ron  tin-  inverse  id! 

-1  - 1 /k  * /° 

x = V (u)  a t b [(l-u)  - 1 | 

where  u is  a uni  tot m (0,1)  variate. 


c.i  nkkai.i:  : i.amupa  nis;::ti.rru’N 

Tukev  l W ] eons  idet a d t ho  inverse  edt 

x = K \u)  ( u - (1-u)  ]/A  (0  ' u ' 1 ) 
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Sclimeis.  t and  Deut.sc 

'll  13'.  | llis.-u: 

;s 

x = K-1(U>  - q 

+ A,  I b(A3. 

A . ) J <0  < 

M - 

«i  < 

r a3 

where  g(\  , A ) 

i - .1 ) 

‘ X3 

if  u *' 

X. 

l (» - V 

if  u > 

A, 

•I 

with  pat ameter  \ ' . distribution  . 

trie  about  .■<  n>  anti  lias  a wide  ran.e  (l*  *’,•)  o! 
values  ot  :•  , . Jo  i net  and  Rosenblatt  | .'  > | noted 
that  lor  A * . 1 tin-  distribution  elosely  approxi- 

mated the  normal  d i.-  tr  ibut  ion. 

Kamb  t rg  and  Selim  iset  ] 12 , 11]  generalized  the  il  i s- 

t r i but  i i a to 

1 \ \ 

x K (w)  \ +[u  -(1-u)  4]/X  (0  u •_  1) 

where  V,,  A , A have  tin-  same  sign.  The  loeat ion 
par. i i ter  is  \^,  the  seale  parameter  is  A,,  and 
shape  is  determined  by  ‘ ^ and  \ . This  distribu- 
tion eovers  all  ot  the  ("j.  ” ,)  plane  « Xt«pt  that 
rot'll  .pondiup.  tt'  tin-  V - shaped  beta  d i s t r i but  ion . 
Comb  i tied  with  the  !’  shaped  beta,  the  entire  l j , 
{',)  plane  mav  in-  eovetvd  in  two  functional  terms 
by  not  i tty,  that  the  limit  in,;  distribution  ot  the 
gemtuli.*ed  lambda  distribution  as  either  \ ..  or  X. 
goes  to  infinity  i ; t!ie  power  futtei  ion  d is  t r i but  ion, 
wbi i b is  a rpoeial  ease  ot  the  beta  family  for  one 
parameter  equal  to  one. 

Nearly  espmit'nt  i a 1 and  nearly  normal  variates  may 
I • enerated  from  th<  generalized  lambda  distrlbu- 
tion.  The  exponential  is  a limiting  distribution 
as  \ | ; \ , L • A , • goes  to  infinity  and  \ 

)’.'’i  to  zero.  The  ?.  t and.  trd  i red  normal  distribution 
is  approx i:. at ed  bv  the  patameter  values  A 0, 

A,  . I / > and  A ^ V .11.0,  wbi  eb  mat  eh  the 
first  tout  : o:  ell  ( -.  flu  maximum  dtviation  between 
the  two  eilt'  ; is  about  .001.  (This  is  about  one- 
tenth  the  maximum  rt  lor  o.btaiueil  using  the  central 
Ii;  it  theoi.  i appro  . i mat  i en  with  twelve  uniform 
variates.)  While  not  ta.t  tAtkin-ni  and  l’eatee  |2]), 
Ibis  doe-,  vield  .i  i»ue  line  normal  getter. itot  lor 
situatior.  where  ah  an  approxinat ion  is  reason- 
able. (lhefe  are  ; • tin  taster  and  ex.ict  t eehn  i ques 
lor  sots  al  .'■•iii  i .it  i.'ii,  whieh  v.*u  probablv  already 
know  it  von  have  re. id  this  far.) 


where  • * X < cu,  0 < A,,  0 < X ^ , and  0 • 1 

whieh  is  a mixture  ot  two  distributions.  To  the 
author's  knowledge,  this  is  i lie  on  I v distribution 
having  a single  funetional  form  whieh  severs  the 
int  ire  (e  , ,)  plane.  Three  d i st  r i but  ions  are  ob- 

tained  exactly:  llernoulli  trials  if  X 0,  uni- 

form if  A I,  and  exponential  as  A go.  . to  zero 

with  A,  A X / \ and  X , 1 . 

J -3  s 

Parameters  mav  be  determined  using  on  1 v elosed  form 
expressions  bv  matching  pereent  i Les.  I.’ i von  the 
mode,  Prob  {x  • mode],  and  any  two  pereentiles  u, , 
Pj)  and  (x,,  1',),  the  correspond  i ng  pa  i . .tit  ate 
Xj  mode,*  X.  ‘ Droit  {x  < mode}, 

1.1|01-X1)/0|-K,)| 

h I II I (.  \/(-|i|  )/  ( \ ,(-p  I 

and  ^ 

V,  ■ |Xj  - xg/p.  - Pjl  3 

Moments  may  be  mat  eh. si  using  graphs  in  Sehmeis.-r  and 
Dent  sell  [ 3-»  | . Random  numbet  generat  ion  is  as. tin 
st  ra  iglit  1 orward  due  to  the  closed  form  invei-e  edt. 

An  important  disadvantage  ol  the  absolute  lambda 
distribution  is  it.-,  t ruinated  tails  and  extreme  va- 
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V.  I c >Ni  I t '- U'N 

Available  l llu-  simulation  practitioner  are  ni.my 
technique  . which  .iIIkw  1 1 it*  mode  l 1 i n • ami  genora- 
t ion  »!  j : -'bah  i ! i si  i . a.  poju-m  i.i  :•  -i  t bo  -ia»dard 
distribution  arr  not  adequate.  Tin*  techniques 
vary  in  enoi  . • . ■ < c»  I pa r nine t e i * t i it  i 

rase  ot  variate  a iiri.ii  i-'ii,  .md  -.pil'd  of  execut  ion . 
Tin*  technique  :■<  lei'toil  should  ret  led  tin  particu- 
lar application  at  hand. 
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